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Abstract 

We prove the integrability of the discretization of the Neumann system re- 
cently proposed by V. Adler. 
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1 Introduction 



The famous Neumann system belongs to the area of constrained mechanics. It de- 
scribes the motion of a particle on the surface of the sphere S = {x G : (x,x) = 1} 
under the influence of the harmonic potential |(fix, x), where Q = diag(a>i, . . . , u>n). 
(Here and below (•, •) stands for the standard Euclidean scalar product in R N ). The 
equations of motion read: 

Xk = -u k x k - ax k , 1 < k < N , (1) 

where a = a(x,x) is the Lagrange multiplier assuring the validity of the relations 
(x,x) = 1, (x,x) = during the evolution. It is easy to see that 

a = (x,x) — (Qx,x) . (2) 

This system can be given a Hamiltonian interpretation in terms of the Dirac Poisson 
bracket, and it turns out to be completely integrable in the Liouville-Arnold sense, 
all integrals being quadratic in momenta (which also implies that it is solvable via the 
separation of variables method). The literature on the Neumann system includes the 
original |J, as well as modern presentations, e.g. @, ||, |[|. 

The present paper is devoted to a discretization of ([[]), introduced recently by 
V. Adler in a beautiful paper ||, along with a novel discretization of the Landau- 
Lifschitz system. In difference equations for x : hZ i— > S, we write x for x(nh), x 
for x((n + l)h), and x for x((n — l)h). The equations of motion of the V.Adler's 
discretization read: 



x k + x k x k + x k (2 - -^-uj k )x k 



~k i -ZK = 2 1< k< N (3) 

l + (x,x) l + (x,x) l-£(nx,x)' - - KJ 

While the integrability of this system for = 3 follows from Adler's results, the 
general case was left open by him. We construct explicitly the full set of integrals of 
motion for ([|). 

It has to be said that several different integrable discretization of the Neumann 
system are currently known. The first one of them, due to A. Veselov 0, |8| is governed 
by the following difference equations: 

x k + x k = (1 + h 2 uj k )- 1/2 f3x k , 1 < k < N , (4) 

where (3 is the Lagrange multiplier, assuring that x remains on the sphere S during 
the discrete time evolution. It is easy to see that 

2((/ + /i 2 ft)- 1 /2 X) x) 2((J + h 2 n)- 1 / 2 x,x) ((I + h 2 ti)- 1 / 2 x } x+ x) 

\i + h 2 ny i x,x) ~ ((i + h 2 ny i x,x) ~ ((/ + h 2 ny i x,x) 



The first two expressions for (3 make this difference equation explicit for the evolution 
in both the positive and the negative directions of (discrete) time. The third expression 
makes it obvious that (3 = 2 — h 2 a + o{h 2 ) with a from (g), so that we indeed 
have a discretization of (|l]). This equation is a discrete time Lagrangian system, 
which allows one to introduce the canonically conjugate momenta in a systematic 



way. Performing this, the main feature of (^) comes into the light: namely, it shares 
the integrals of motion with its continuous time counterpart. In other words, it is a 
Backhand transformation for the continuous time Neumann system (cf. [§]). 

The second integrable discretization of the Neumann system, discovered by O. 
Ragnisco [1C], is governed by the following difference equations: 

2xk + ~ fc = —h 2 uJkXk + h 2 (Qx,x)xk , 1 < k < N . (6) 



(x, x) (x, x) 

Again, this is a discrete time Lagrangian system. Introducing the canonically con- 
jugate momenta in a proper way, one sees that the integrals of motion for (^) are 
0(/i)-perturbations of the integrals of the continuous time Neumann system. 

The comparison of the previously known discretizations was performed in |l T| . 



With the advent of the V. Adler's discretization and the proof of its integrability, 
the Neumann system becomes the best candidate for the championship in possessing 
various different integrable discretizations. 

2 Neumann system 

The Neumann system describes the motion of a point x G R N under the potential 
h(Qx, x), constrained to the sphere 

S = {xeR N : (x,x) = l}. (1) 

The Lagrangian approach to this problem is as follows. The motions should deliver 
local extrema to the action functional 



ti 

to 



L(x(t),x(t))dt 



where L : TS i— > M. is the Lagrange function, given by 

L(x,x) = ~(x,x) - ^(Qx,x) - ^a((x,x) - 1) . (2) 

Here the first two terms on the right-hand side represent the unconstrained Lagrange 
function, and the Lagrange multiplier a has to be chosen to assure that the solution 
of the variational problem lies on the constrained manifold TS, which is described by 
the equations 

(x, x) = 1 , (x,x) = . (3) 
The differential equations of the extremals of the above problem read: 

X k = -UkXk - OLX k , (4) 

or, in the vector form, 

x = —Qx — ax . (5) 
The value of a is determined as follows: 

= (x,x)' = (x, x) + (x, x) = (x, x) — (Qx, x) — a . 



Therefore 

a = (x,x) — (ftx,x) . (6) 

So, the complete description of the Neumann problem is delivered by the equations 
of motion @, augmented with the expression (||). 

The Legendre transformation, leading to the Hamiltonian interpretation of the 
above system, is given by: 

K(x,p) = (x,p) — L(x, x) , 

where the canonically conjugate momenta p are given by 

p = dL/dx — x . (7) 

Hence we obtain ^ ^ 

H(x,p) = -x{p,p) + ^x,x) . (8) 

The corresponding symplectic structure is the restriction of the standard symplectic 
structure of the space R 2N (x,p), 

{Pk,Xj} = 5 kj , (9) 

to the submanifold T*S which is singled out by the relations 

<j> 1 = (x,x)-l = 0, <t> 2 = <p, x) = . (10) 

The Dirac Poisson bracket for this symplectic structure on T*S is given by the fol- 
lowing relations: 

{x k ,Xj} D = 0, {pk,Xj}u = hj ~ t^t , {Pk,Pj}D = ^-j — t^- ■ (11) 

{x,x) {x,x) 

It is easy to check that the Hamiltonian vector field generated by the Hamilton func- 
tion (d) in the bracket {-, -}d on T*S, is given by: 

x k = {H.,x k } D = p k , p k = {R,p k } D = -u} k x k - ax k , (12) 

where 

a = (p,p) — (Qx,x) . (13) 

Obviously, this is nothing but the first-order form of (|4]) with the multiplier (0). 

Supposing that all u>k are distinct (which will be assumed from now on), one can 
prove that the following N functions are integrals of motion of (|P2"D: 

F fc = xl + £ {PkXj " XkPj)2 , l<k<N. (14) 

j^k U k - UJj 

Only N — 1 of them are functionally independent on T*S, due to the relation 

N 

X)F fc =(x,x), (15) 

k=l 

and this is equal to 1 on T*S. The Hamilton function of the Neumann system may 
be represented as 

1^1 1 
H = -J2^ k F k = - (jp,p)(x,x) - (p,x) 2 ) + -(ttx,x) , (16) 

which coincides with on T*S. All functions are in involution, assuring the com- 
plete integrability of the Neumann system with respect to the Dirac Poisson bracket 
{v}d. 



I 



3 Adler's discretization 



h 2 

-0. 

4 



(17) 



It will be convenient to denote 

h 2 

a k = 1 - — u k , A = diag(«i, . . . , a N ) = I 

In these notations, the Adler's discretization (|3]) is written as 

x k -\- x k ^ x k -\- x k 2cY k x k (18) 

1 + (x, x) 1 + (x, x) (Ax, x) 

From the first glance, these equations look implicit for the evolution in both the 
positive and the negative directions of time. However, this is actually not the case. 
Consider, for definiteness, the case of the positive time direction. We can express 
from (|i~8D the quantity 1 + (x, x) through x and x. Indeed, rewrite ([L8]) as 

x + x 2Ax x + x 

l + (x, x) (Ax,x) l + (x,x) 

Taking the square of the norm of both sides, we find: 



1 + (x,x) 



2Ax 



X + X 



(Ax, x) 1 + (x, x) 

With this at hand, we can use (|i"9|) to express also x as function of x and x. 

By the way, the previous argument has an important corollary. It is easy to see 
that the last formula may be represented as 

2(Ax,Ax) 2(Ax,x) 



1 



1 



+ 



l + (x,x) 1 + (x, x) (Ax,x) 2 (Ax, x) (1 + (x, x)) 

However, we could as well perform a similar calculation for the evolution in the neg- 
ative direction of time, which would lead to the following formula: 



1 



1 



+ 



2 (Ax, Ax) 



2 {Ax, x) 



(20) 



(21) 



1 + (x, x) 1 + (x, x) (Ax,x) 2 (Ax, x) (1 + (x, x)) 

Comparing the last two formulas, we come to the following conclusion: 

(Ax, x) (Ax, x) 
1 + (x, x) 1 + (x, x) 

In other words, we found an integral of motion of the difference equation flT8"|): 

2(Ax,x) 
h{x,x) = . 

1 + (x, X) 

The second integral of motion can be obtained even more straightforwardly. Indeed, 
from (|i~9D there follows easily: 

(A _1 (5 + x), x + x) _ (A~ x (x + x), x + x) 

(l + (x,x)) 2 ~ (l + (x, x)) 2 ' ( j 

In other words, the following function is also an integral of motion of the difference 
equation (Jl8|): 

(A _1 (x + x), x + x) 



l2[X,X) 



[l + (x,x)Y 



(23) 



c 



4 Lagrangian formulation 

We now demonstrate that ([18]) may be interpreted as the equation of extremals of a 
discrete time action functional 

m 

s= £ i>(x((n + l)h),x(nhn , 

n=no 

where LiSxSwRisa discrete time Lagrange function. Recall that the discrete 
time Lagrangian equations of motion read: 

V x h(x, x) + V x h(x, x) = , (24) 

and that the momenta p G T*S, p G T~S canonically conjugate to x, resp. to x, are 
given by: 

hp = -V x h(x,x) , (25) 

hp = V~L(x, x) . (26) 

Consider the discrete time Lagrange function on S x S: 

hL(x,x) = -2 log(l + (x,x)) + 2 log(Ax,x) , x G S 1 , x & S . (27) 



Applying (p5|), fl26|), we find: 



2xk 4=OikXk 
1 + {x, x) {Ax, x) 

fcfc = -n^-v+^t- ( 29 ) 

1 + (x, x) 

Here the scalar multipliers 7, 5 have to be chosen so as to assure that 

P et:s, peT^S, 

or, in other words, to assure that the following relations hold: 

(p,x> = 0, (p,x)=0. (30) 
It is easy to see that this is achieved if 

7= ^ v -4= L-,-2, f = 2( ^> = ^ + 2. (31) 

l + (x,x) l + (x,x) 1 + (x, x) 1 + (x, x) 

So, the following are the Lagrangian equations of motion of the Adler's discretization 
of the Neumann system: 

hpk = ; * fc ; ---^r + 2x fc , 32 

1 + (x, x) (Ax, x) 

, ~ 2(xi. + Xfc) _ . , 

hp k = - /* \ +2£ fc . 33 
1 + (x, x) 



Clearly, these two equations yield also the Newtonian form ( |T5| ) of the equations of 
motion. 

The equations fl32[) , C|33|) define a symplectic map (x,p) G T*S i— > (x,p) G T*,S 
explicitly. Indeed, the first equation fl32|) yields 



2(x fc + x fc ) 



4a k x k 



which implies 



1 + (x, x) (Ax, x) 
2 



- 2x k + hp k 



2Ax hp 
— x ~\~ — 

(Ax, x) 2 



(34) 



1 + (x, x) 

This, substituted back into (j32p, allows us to determine x, and then, finally, the second 
equation of motion ( |33"D defines p. (Clearly, this is a perifrase of the similar argument 
in the previous section). 



5 Integrability 

We can now express the integrals of motion Ii, I2 in terms of canonically conjugate 
variables (x,p). Straightforward calculations verify the following formulas: 

h(x,x) = 2 ( Aa 2 x ) = (Ax, a;) - h(Ax,p) - ^-(p,p)(Ax, x) , (35) 
1 + (x, x) 4 

/ 2 (x,x) = (A ' 1 { 3 ; + / ! ;), ^ 2 +X) = (A-^x) - h(K-\,p) + ^(A-^p) . (36) 
(1 + (x,x)J^ 4 

These two integrals are enough to assure the Liouville- Arnold integrability of the 
map (^3), (^) for N = 3. The full set of integrals in the general case is given in the 
following statement, which constitutes the main result of this Letter. 

Theorem 1 // all a k 's are distinct, then the following functions are integrals of mo- 
tion of the Adler's discrete time Neumann system: 

r 2 u . h2 ( x kPj ~ Xjp k )(a k x k pj - ajXjPk) 

T k = x k - hx k p k + — > J - J - J - . (37) 

4 ^ atj - a k 

Proof. Denote 

X fc j = x k pj - xjp k , Y kj = a k x k pj - ajXjPk . (38) 

We have: 

Using the equations of motion fl32"D , (p3|), we find: 



2 J 1 + (x, x) 



^ , , a k x k Xj - ajX k Xj + (a k - aj)x k Xj , 

2 1 + (x, x) 



and 



h — 2xkXj 



2 aij — a k (Ax, x) 

kj ^ kj Xj^Xj ~\~ XfcXj ~\~ XfcXj ~\~ Xj~Xj 

2 1 + (x, x) 
Calculating further, we find: 

4 j£ a, } -a k 

1 ^ — ^ , s. f ^k^j ~ I - %k^j H - %k%j H - ^k^j 



H~~ ^k^j ^k^j • 



E, ~ ~ \( %k%j H~ %k%j H - %k%j H - %k%j ~ ~ \ 

6*^(1 + (x, sc)) - 2x /c x/ c — 5^(1 — (5, x 



1 + (x, a?) 

X k X k h%kPk j 



4 j^fe a i ~~ a fe 



X}-Xj I i /~ \ CHj jXfcXj 



(Ax, x) V l + (x,x) 

2x fc / (x fc + x fc ) \ 

= tt r 2a fc x fc - - — - (Ax, x) - x k (Ax, x) 

{Ax, x) \ 1 + {x, x) J 

= -hx k pk ■ 
Collecting all results, we find: 

This proves the Theorem. I 

The involutivity of the integrals may be established by the arguments similar 
to those in [Q. It is easy to see that the two previously found integrals 
admit the following expression in terms of T k . 

N N 

Hi = a k^k , h = a k Fk ■ 

k=l k=l 



6 Conclusions 

Several things remain to be done to put the V. Adler's discrete time Neumann system 
into the modern framework of discrete integrable systems: to find Lax representations 
and their r-matrix interpretation, and to perform an integration in terms of theta- 
functions. I hope to have an opprtunity to present these results elsewhere. I thank 
V. Adler and A. Shabat for showing me their results prior to publication. 
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